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Lattice calculations of the form factors for the charm semileptonic decays D — > Klv and D — > nlv 
provide inputs to direct determinations of the CKM matrix elements \V CS \ and \V a /\ and can be 
designed to validate calculations of the form factors for the bottom semileptonic decays B — > nlv 
and B — > Kll. We are using Fermilab charm (bottom) quarks and asqtad staggered light quarks on 
the 2+1 flavor asqtad MILC ensembles to calculate the charm (bottom) form factors. We outline 
improvements to the previous calculation of the charm form factors and detail our progress. We 
expect our current round of data production to allow us to reduce the theoretical uncertainties in 
| V -| and \V cd \ from 10.5% and 11%, respectively, to about 7%. 
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1. Introduction 

The CKM matrix elements |V CJ | and |V c d| can be extracted to greatest precision (currently to 
0.02% and 0.4%, respectively) by assuming CKM unitarity and performing a fit to all data M\. 
However, the simplest tests of unitarity require direct determinations of the CKM matrix elements. 

The decay rate for D — > K(%)lV is proportional to a form factor and |V„| (|V C </|). Experiments 
can measure the decay rates and the form factor shapes, but nonperturbative calculations of the 
strong force are required to fix the form factor normalizations and extract \V CS ^) I- Therefore these 
decays allow direct determinations of \V cs u) | and consistency checks between lattice QCD and 
unitarity. Such consistency increases our confidence in both. 

In June CLEO-c published the results of an analysis of 818 pb _1 collected at charm threshold 
[g]. Combining the CLEO-c results with the first 2+1 flavor lattice calculations of the D — > K(n)lv 
form factors |] yields \V cs(d) | |@|: 

|V„| = 0.985(1 ±0.9% ±0.6% ±10.5%), (1.1) 
\V cd \ = 0.234(1±3%±0.9%±U%). (1.2) 

The first errors are experimental statistical errors, and the second are experimental systematics. The 
third errors are due to uncertainties in the lattice QCD calculations. The theory errors dominate the 
uncertainties. 

Discretization effects are the dominant source of the theory errors [||]. Other uncertainties 
enter because of incomplete suppression of oscillations due to opposite-parity states, truncation 
effects in fits to staggered chiral perturbation theory (S#PT), and model-dependence implicit in the 
Becirevic-Kaidalov (BK) parameterization 

These sources of uncertainty were addressed in work on B — > nlv decays By calculating 
the D — > K{n)lv form factors using the same methods, we may be able to validate their application 
to calculations of the form factors for B — > nlv and B — > KlT. The former decay allows a precise 
determination of \V u b\ and a stringent test of unitarity. The latter is a rare decay and a prime 
candidate for new physics. Below we describe our progress in reducing the uncertainties in the 
charm form factors and anticipate the reduction of the uncertainties in | V cs /^ \. 

2. Ensembles and quark masses 

To decrease discretization effects and improve our control of the chiral extrapolation, we are 
generating full QCD and partially quenched data on each of the ensembles shown in Table [j]. 
These ensembles include the four most chiral coarse ensembles used in the calculations of Ref. 
the two fine (a as 0.09 fm) ensembles included in our recent calculation of the form factor for 
B — > nlv two additional fine ensembles, three superfine (a w 0.06 fm) ensembles, and one 
ultrafine (a « 0.045 fm) ensemble [0]. The MILC Collaboration has increased the number of 
configurations in each of the previously used coarse and fine ensembles by a factor of four, and we 
expect a corresponding decrease in all statistics-dominated uncertainties by a factor of two. 

We have found that randomizing the spatial location of the sources significantly decreases 
autocorrelations in 2-point functions, which suggests that we may be able to increase our statistics 
further by increasing the number of source times on each configuration. We have nearly completed 
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a (fm) 


am\ 1 am s 


Volume 


Neon f 


cini va \ ence 


coarse 


0.12 


0.02/0.05 


20 3 x 64 


2052 


0.005, 0.007, 0.01, 






0.01/0.05 


20 3 x 64 


2259 


0.02, 0.03, 0.0415, 






0.007/0.05 


20 3 x 64 


2110 


0.05; 0.0349 






0.005/0.05 


24 3 x 64 


2099 




fine 


0.09 


0.0124/0.031 


28 3 x 96 


1996 


0.0031, 0.0047, 0.0062, 






0.0062/0.031 


28 3 x 96 


1946 


0.0093, 0.0124, 0.031; 






0.00465/0.031 


32 3 x 96 


983 


0.0261 






0.0031/0.031 


40 3 x 96 


1015 




superfine 


0.06 


0.0036/0.018 


48 3 x 144 


668 


0.0036, 0.0072, 0.0018, 






0.0025/0.018 


56 3 x 144 


800 


0.0025, 0.0054, 0.0160; 






0.0018/0.018 


64 3 x 144 


826 


0.0188 


ultrafine 


0.045 


0.0028/0.014 


64 3 x 192 


861 


TBD 



Table 1: Asqtad staggered quark ensembles generated by the MILC Collaboration [| ^J] and slated for 
upcoming heavy-light analyses, together with the valence quark masses being used at each lattice spacing. 
The last valence mass listed at each lattice spacing (after the semicolon) is the tuned strange quark mass. We 
are presently generating correlators at four source times on each ensemble and investigating the possibility 
of adding more source times to further increase the total number of source-configurations. 

data generation at four source times on the coarse ensembles, the fine ensembles with mi = 0Am s , 
0.2m s , and 0Am s , and the superfine ensemble with mi = 0.2m s . 

Power counting arguments [J|, ^] indicate that including these ensembles will effectively elim- 
inate discretization effects due to light quarks and gluons, while heavy-quark discretization effects 
will be reduced but remain significant. To improve our estimates of heavy-quark discretization 
effects, we are investigating including them in chiral-continuum expansions JTfj|]. This approach 
incorporates the information from power counting while more systematically fixing the appropriate 
hadronic scales. 

3. Correlators and correlator ratios 

The form factors parameterize the hadronic matrix elements of the flavor-changing vector 
currents, 

(K{n)\V^\D) = y/2m^ 

where Vu is the lattice current corresponding to isy^c {idy^c), v = Pd/^id is the four-velocity of 
the D meson, p± = p K ^ — {pk(tl) ' v ) v is ine component of kaon (pion) momentum perpendicular 
to v, and q 2 = (pp — Pk(k)) 2 i s ine invariant mass of the leptons. We work in the D-meson rest 
frame, in which the form factors are proportional to the temporal and spatial components of the 
hadronic matrix elements, and q 2 = m 2 D + m 2 K ^ — 2moE K ( K y 

One way to extract the hadronic matrix elements is by considering simple ratios of 3-point to 
2-point correlators [||], 

C^(t,T; PK{7[) ) 
C^\t; Vm )C°{T-t) 



r 1 



I 



(q )+p±nf ± 



(r 



(3.1) 
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where T is the separation between source and sink in the 3-point functions, and 

c^ K{K \t,T-p m ) = £^w y (^ w (^,o)y M (f, y )4( ?/ , x )), 

? G [*;, f/ = (ti + r) mod n t ], (3.3) 
cf W (np* W ) = I^ w ' x (^ W (^0)^ (7t) (r,x)), 

X 

t G [f,-, f/ = {ti + n,) mod n t ), (3.4) 
C? (0 = £<^dM)€^,x)), t G [?,-, f/ = {ti + n,) modn,). (3.5) 

X 

where «, is the temporal extent of the lattice, and is the momentum of the outgoing kaon 
(pion). We calculate the correlators for momenta Pk(tl) = (0,0,0), (1,0,0), (1, 1,0), (1,1, 1), and 
(2,0,0) (in units of 2n/L, where L is the spatial extent of the lattice) and all times t in the ranges 
shown. We increase statistics by averaging correlators with source times f; = 0, n t /4, n t /2, 3n t /4. 
The D-meson interpolating operators G are smeared with a charmonium wavefunction to suppress 
coupling to excited states. 

C3 is calculated with insertions of the current operator at all times t between the source and 
sink. At sufficiently large source-sink separations T and times t sufficiently far from both source 



and sink (0 <C t <C T), a plateau emerges in the ratio (3.2). This plateau is directly proportional to 
the desired hadronic matrix element. 

In practice we find that oscillations from opposite-parity excited states contaminate the en- 
tire plateau region ^ |6|]. We therefore consider the more carefully constructed correlator ratios 
introduced in Ref . [Hi : 



C?~'(t;p m )(Z(T-t) 

where = (1; Pk(x)) an d tne correlators C3, C2 are constructed from the correlators C3, C2 

to eliminate oscillations from opposite-parity states: 



C 3 (t, T ) = - |c 3 (t, T) + C 3 (t , T + 1 )e m ° + 2C 3 (t + 1 , F)Aw ~ m ° + 2C 3 (f + 1 , T + 1)A 
8 L 

+ C 3 (f + 2, r)e 2 ^W ~'" d) + C 3 (f + 2, T + 1 )e 2£ *w ~ m ° 



(3.7) 



C 2 (t) = i [C 2 (0 +2C 2 (f + l)e mD + C 2 (f + 2> 2mD ] . (3.8) 
Experience suggests that the errors in direct fits to the oscillating states can be larger than errors in 



simpler fits. The construction of (3.6) and (3.7, 3.8) allows us to fit the ratios to constants without 



introducing systematic errors. In the plateau region (0 <C t <C T), the ratios R%7^ for ji. = 

(jU = approach the form factors (f±^ K ^)- 

For source-sink separations T = 16 and T = 20, examples of the plateaus are shown in Figs. |T] 
and Q where the features leading to the choice of these T-values can also be seen. As the source- 
sink separation increases, signal-to-noise decreases. As the source-sink separation decreases, the 
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Figure 1: Ratios of correlators for extracting the form factor f®^ K {q 2 )- The correlators were calculated on 
2110 configurations of the coarse ensemble with m; = 0.14mj. T — 16, 20 are the source-sink separations, 
and the three-momenta p of the pions are given in units of 2n/L, where L is the spatial extent of the lattice. 
Note the excited state contamination in the zero momentum data with T = 16. 

plateau region shrinks and eventually disappears. The optimal T-value is the smallest for which a 
plateau exists. For this T, signal-to-noise is maximized without sacrificing the plateau to excited 
state contamination. The statistical errors increase with momentum, so the optimal T is momentum 
dependent. 

To optimize T we generated data with T = 16, 18, and 20 on the coarse m/ = 0.14m s ensem- 
ble. As shown in Figs. [I] and for T = 20 plateaus exist for all momenta. At zero momentum, 
comparing the T = 16 data with the T = 20 data reveals the effects of excited state contamination 
in the T = 16 data for all t; the plateau has essentially vanished. At nonzero momentum, comparing 
the T = 16 data with the T = 20 data reveals smaller statistical errors in the T = 16 data with intact 
plateau regions. The larger T allows checks for excited state contamination at smaller momenta, 
and the smaller T allows us to minimize statistical errors at larger momenta. On the remaining en- 
sembles, we expect the optimal T-value in physical units to be similar. We are therefore generating 
data on each ensemble with two aJ-values of approximately 0.12 fm x 16 and 0.12 fm x 20. 

4. Renormalization and chiral-continuum-energy extrapolation-interpolation 

Lattice form factors obtained from the plateaus in Figs. [I] and ||] must be renormalized and 
extrapolated to zero lattice spacing and the physical light quark masses. The renormalization factors 
can be written as products of non-perturbatively calculable factors Zy and perturbatively calculable 
factors p. The uncertainties in these renormalization factors contribute to the uncertainties in the 
form factors and CKM matrix elements. 

To perform simultaneous chiral-continuum extrapolations and the kaon (pion) energy inter- 
polation, we can use staggered heavy meson partially quenched chiral perturbation theory (#PT) 
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Figure 2: Ratios of correlators for extracting the form factor ff~^"(q 2 )- The correlators were calculated on 
2110 configurations of the coarse ensemble with to/ = 0.14m s . T — 16, 20 are the source-sink separations, 
and the three-momenta p of the pions are given in units of 2%/L, where L is the spatial extent of the lattice. 
The consistency of the results for T = 16 and T = 20 indicates that the smaller source-sink separation can 
be used to minimize statistical errors without introducing significant excited state contamination. 



with constrained curve fitting [jS, 11, [12]. This approach incorporates the energy-dependence of the 
form factors and yields a model-independent result while accounting for the systematic error due 
to truncating the expansion. 

To extract [V^^ |, one can divide the experimental results [^] by the lattice form factors eval- 
uated at q 2 = 0. However, minimizing the uncertainty in \V cs u) j requires a simultaneous fit to all 
(experimental and lattice) data. The analyticity-based parameterization described in Ref. [13] cap- 
tures the energy-dependence of the form factors throughout the kinematic domains, so using it to fit 
the data and extract CKM matrix elements does not introduce model-dependent systematic errors. 

For D — > K(n)lv, the energy-domains of the lattice and experimental data overlap signifi- 
cantly, allowing a stringent test of the consistency of the shapes of the form factors as determined 
independently by the lattice and experiment. This test provides important validation for applying 
the analyticity-based parameterization to the extraction of \V u b \ from B — > nlv, in which the overlap 
of the lattice data and experimental data is smaller and this self-consistency check, less powerful. 



5. Expected uncertainties 

A projected error budget for the form factors at q 2 = is shown in Table ||. The expected 
uncertainties reflect previous experience with B — > nlv including the use of improved correlator 
ratios, #PT with constrained curve fitting, and the analyticity-based parameterization to eliminate 
systematic errors due to incomplete cancellation of oscillating state contributions, truncation of the 
chiral expansion, and model-dependence in the BK parameterization. The projections also reflect 
the four-fold increase in statistics on the coarse ensembles and the addition of the two largest 
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Stat. + ^PT 


§d*dtc n 


m 




HQ 


Zy 


p L 3 < oo 


Sys. 


Total 


4.9 


2.9 1.4 


0.3 


1.3 0.2 0.1 


3.9 


0.7 


0.7 0.5 


5.4 


7.3 



Table 2: Contributions to the relative uncertainties in the form factors at q 2 — assuming data with four 
source times on the four extended coarse ensembles, two largest fine ensembles, and the superfine m/ = 
0.2m s ensemble. The errors are due to limited statistics and the truncation of chiral perturbation theory; 
uncertainties in the D*Dn coupling, scale, average up-down quark mass, strange quark mass, and charm 
hopping parameter; momentum-dependent discretization effects from the light quarks and gluons; heavy- 
quark discretization effects; uncertainties in the renormalization factors Zy and p; and finite volume effects. 
The last two entries are the total systematics and the total error, both added in quadrature. 

fine ensembles and the superfine mi = 0.2m s ensemble. The increase in statistics decreases our 
statistical uncertainties by a factor of two, while the addition of the superfine ensemble reduces 
systematic errors due to heavy-quark discretization effects. 

Heavy-quark discretization effects and the uncertainty in the D*Dn coupling dominate the 
systematic uncertainties, while statistics and #PT truncation error are alone comparable to the 
entire remaining systematic error. Heavy-quark discretization effects are sensitive to the smallest 
lattice spacings included, so they will decrease further with the addition of the ultrafine ensemble 
m Tabled The error due to the D*Dn coupling may respond to the increased statistics. From 
Table ^| and Eqs. Ql.l[ ) and (|L2|), we expect to reduce the theoretical uncertainties in the CKM 
matrix elements from about 11% to about 7%. 

Fermilab is operated by Fermi Research Alliance, LLC, under Contract No. 
DE-AC02-07CH11359 with the United States Department of Energy. 

References 

[1] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667, 1 (2008). 

[2] D. Besson et al. [CLEO Collaboration], Phys. Rev. D 80, 032005 (2009) [arXiv:0906.2983 [hep-ex]]. 

[3] C. Aubin et al. [Fermilab Lattice and MILC Collaborations], Phys. Rev. Lett. 94, 01 1601 (2005) 
[arXiv:hep-ph/0408306]. 

[4] C. Bernard et al. [Fermilab Lattice and MILC Collaborations], Phys. Rev. D 80, 034026 (2009) 
[arXiv:0906.2498 [hep-lat]]. 

[5] D. Becirevic and A. B. Kaidalov, Phys. Lett. B 478, 417 (2000) [arXiv:hep-ph/9904490]. 

[6] J. A. Bailey et al. [Fermilab Lattice and MILC Collaborations], Phys. Rev. D 79, 054507 (2009) 
[arXiv:08 11.3640 [hep-lat]]. 

[7] A. Bazavov et al, arXiv:0903.3598 [hep-lat]. 

[8] C. W. Bernard et al. [MILC Collaboration], Phys. Rev. D 64, 054506 (2001) [arXiv:hep-lat/0 104002]. 

[9] C. Aubin et al. [MILC Collaboration], Phys. Rev. D 70, 094505 (2004) [arXiv:hep-lat/0402030]. 

[10] A. Bazavov et al. [Fermilab Lattice and MILC Collaborations], PoS(LAT2009)249. 

[11] C. Aubin and C. Bernard, Phys. Rev. D 73, 014515 (2006) [arXiv:hep-lat/0510088]. 

[12] C. Aubin and C. Bernard, Phys. Rev. D 76, 014002 (2007) [arXiv:0704.0795 [hep-lat]]. 

[13] T. Becher and R. J. Hill, Phys. Lett. B 633, 61 (2006) [arXiv:hep-ph/0509090]. 



7 



